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(57) Abstract 

Provided is a method of measuring impedance based on carrier function Laplace transform, comprising the steps of: (1) detecting a 
response signal from an object device to which signal excitation such as pulse, interrupt or constant load is applied; (2) obtaining parameters 
of a carrier function by linearly or non-lineariy fitting the response signal in step (1) to the carrier function allowing an analytical Laplace 
transform; (3) calculating an impedance function in the Laplace domain by using an analytic relation between the carrier function and the 
Laplace transfer function calculated from carrier function obtained by the fitting results in step (2); (4) calculating the impedance spectrum 
in a frequency domain by using the parameters of the Laplace domain impedance function calculated in step (3), and calculating the 
measurement error of the frequency domain impedance spectrum using the standard deviation of the parameters and its correlation; and (5) 
displaying the frequency domain impedance spectrum calculated in step (4) and storing the result. 



BNSDOCID: <WO 0055639A1_I_> 



WO 00/55639 



PCT/KR99/00815 



LAPLACE TRANSFORM IMPEDANCE SPECTROMETER AND ITS 

MEASUREMENT METHOD 



5 Technical Field 

The present invention relates to a Laplace transform impedance spectrometer 
and its measurement method. More particularly, the method of and apparatus for 
measuring Laplace transform impedance can provide a high quality of wide-ranged 
impedance spectrum measurement within a minimal time with, applicable to the 
10 various electrical circuits, non-linear devices, capacitors, and other electrochemical 
devices such as primary and secondary batteries and fuel cells. 

Background Art 

The technique of measuring impedance spectrum of the electrical or* 
15 electrochemical device can be widely applied to a characterization of an electrical 
circuit, material evaluation, corrosion protection, evaluation of properties, quality 
control and capacity estimation of batteries. 

As a widely used measurement for impedance spectrum, there is a method of 
applying a periodic perturbation signal at a certain frequency, and measuring the 
20 amplitude and phase of the response, for example, by using frequency response 
analyzer. 

U.S. Patent Nos. 4,196,475 and 3, 634, 760 introduce a method of analyzing a 
response characteristic at a single frequency. But the method stated in the above, 
patent applications has problems that it must use expensive devices such as a signal 
25 generator and phase detector, and it needs at least two periods of signal to remove 
transient effects. In addition, it takes a long time to proceed the successive 
measurements for each frequency when a spectrum at a plurality of frequencies is 
required. 
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When using the Fast-Fourier Transform (hereinafter, it is referred to FFT) in 
which a response signal is fast-Fourier transformed by means of an excitation signal 
of the multiple-superposed sine waves (refer to G. S. Popkirov and R. N. Schindler, 
Rev. ScL Instrum. 63, 5336 (1992)), the phase detector is not required and the time 
5 required for the spectrum measurement with respect to the plurality of frequencies is 
relatively short in contrast to the method of using the frequency response analyzer. 

Furthermore, the FFT method examines the linearity of the system for testing 
by comparing input and output power spectrums and thus provides a high quality of 
impedance spectra. The FFT method, however, takes more than double the 
10 measurement time corresponding to the period of the minimum frequency in the 
measurement frequency range and, besides, requires a complicated signal generator 
and a large capacity of memory device. 

Since the applied perturbation signal is generated by superposition of non- 
overlapping frequencies, the impedance spectrum is obtained at odd times of the 
15 minimum frequency. Therefore, only a few number of impedance data points can be 
obtained at low frequency region. 

U.S. Patent No. 5,633,801 describes a method of using current pulse instead 
of a plurality of superposed sine waves in impedance measurement. In this method, 
however, considerable frequency dependent noise due to a frequency aliasing and 
20 transient effect sometimes cause a measurement error, thus having a problem in a 
reliability of the measurement results and confidence interval (W. J. Thomson and J. 
R. Macdonald, Proc. Natl. Acad. ScL USA, 90, 6904 (1993)). 

A method stated in U.S. Patent No. 5,794,008, which uses a discrete Laplace 
transform in analyzing the pulse signal in order to get a high quality of measurement 
25 result at a frequency region, requires a long operation time and a large memory size. 

The above method needs, in analysis, computation at each measurement data 
point, and the number of the measurement data points is determined as double of the 
ratio of the maximum and minimum frequencies so that in case of eight decades of 

2 
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measurement sections, one million times of computation process is required for one 
megabyte memory. 

Integration step required for transformation sometimes gives oscillating 

results. 

It is also difficult to determine the accuracy of the measurement using this method, as 
in the pulse FFT, and if the system for testing is non-linear, it is possible to get the 
impedance spectrum from the measurement but there is no physical meaning on it. 

Disclosure of Invention 

Accordingly, an object of the present invention is to provide a carrier function 
Laplace transform impedance spectrometer and its measurement method, obtaining a 
high quality of wide-ranged impedance spectrum within a minimal time for 
characterization of the various electrical circuits, non-linear devices, capacitors, and 
other electrochemical devices such as primary, secondary batteries and fuel cells. -? 

To accomplish an object of the present invention, a Laplace transform 
impedance spectrometer, includes: 

a galvanostat/potentiostat for applying constant voltage, constant current or 
constant load to the electrochemical device by employing constant voltage source, 
constant current source or constant load element such as a resistor and thereby 
detecting voltage and current of the electrochemical device resulting from above 
excitation; 

a voltage/current output unit for removing noise out of voltage and current detected in 
the galvanostat/potentiostat, also removing bias voltage and current and outputting 
signals; 

a two channel analog/digital converter for converting the voltage and current 
output from the voltage/current output unit into digital signal; and 

a control means for fitting the digital voltage and current output from the 
analog/digital converter to Laplace transform carrier function to thereby generate 

3 
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impedance spectrum and storing its data. 

To accomplish the other object of the invention, a method of measuring 
Laplace transform impedance, includes the steps of: 

(1) detecting a response signal of an object whose impedance will be 
5 measured; 

( 2 ) obtaining parameters of a carrier function by linearly or non-linearly 
fitting the response signal in step (1) to the carrier function which has an analytical 
Laplace transform; 

( 3 ) calculating an impedance function in the frequency domain by using an 
10 analytic relation between the carrier function in time domain and its Laplace 

transform function, using parameters of the carrier function obtained through the 
fitting results in step (2); 

{ 4 ) calculating the impedance spectrum in a frequency range determined by 
measurement interval and sampling rate by using impedance function calculated in 
15 step (3), and calculating the measurement error of the frequency domain impedance 
spectrum out of the standard deviation of the parameters and its correlation; and 

(5) displaying the frequency domain impedance spectrum calculated in step 
(4) in a fixed form and storing the result. 

Specifically, the invention provides a Laplace transform spectrometer 
20 providing high measurement accuracy, speed and a simple architecture. 

The invention relates to a method of obtaining impedance spectra by using 
Laplace transform carrier functions suitable for transient response of step signals. In a 
measurement, the spectrometer requires just a small size memory for computation 
and requires only linear operations on measured data which do not include integration. 
25 And, it is different from other art in that it allows use of a standard statistic means for 
assigning confidence interval to a frequency domain measurement data, providing 
therefore analysis of the quality of the measurement and the linearity of the system 
for testing. 
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Brief Description of Drawings 

FIG. 1 is a block diagram showing the structure of a Laplace transform 
impedance spectrometer of the invention; 
5 FIG. 2 is a signal flowchart showing a method of measuring the Laplace 

transform impedance of the invention; 

FIG. 3 shows an impedance spectrum obtained by Laplace-transform of the 
voltage measured in constant load method; 

FIG. 4 is a graph comparing an impedance spectrum obtained by carrier 
10 function Laplace transform of a response signal of a dummy cell consisting of 
resistors and capacitors measured by a current pulse method with an impedance 
spectrum measured by a fast Fourier transforming method; 

FIG. 5 shows an impedance spectrum obtained by carrier function Laplace- ¥ » 

transform of a response signal of a half cell measured by a current interruption r - 
1 5 method in a charged condition; 

FIGS. 6a and 6c are graphs showing the result of an impedance spectrum 
obtained by carrier function Laplace-transform of a response signal of a primary 
alkaline manganese battery in a full charged condition measured by a current 
interruption method and its error-analysis; and 
20 FIG. 7 is a graph showing an impedance spectrum obtained by carrier 

function Laplace-transform of a response signal from a secondary lithium ion battery 
in a full charged condition measured in a current interruption method. 

Best Mode for Carrying out the Invention * 'rr % 

25 Hereinafter, a Laplace transform impedance spectrometer according to the * 

present invention will be described with reference to the accompanying drawings. 

FIG. 1 shows a Laplace-transform impedance spectrometer of the invention. 
As shown in FIG. 1, the present invention requires an electrical or 

5 
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electrochemical device 10 for testing such as electrical circuits, linear device, non- 
linear device, capacitor, primary and secondary batteries and fuel cells; a 
galvanostat/potentiostat for applying constant-voltage, a constant-current or constant- 
load to the electrical or electrochemical device 10 by employing a constant-load 
5 machine such as a constant-voltage source, constant-current source or resistor and 
thereby detecting the voltage or current variation of the electrical or electrochemical 
device 10, resulting from the application of the constant voltage, current or load; a 
voltage/current output unit 30 for removing noises of the voltage or current detected 
by the galvanostat/potentiostat 20, and also of the bias voltage or current and 

10 outputting the signal; a two-channel analog/digital converter 40 for converting the 
voltage or current output from the voltage/current output unit 30 into a digital signal 
in a desired sample extraction speed; and a control means 50 for fitting the digital 
voltage and current output from the analog/digital converter 40 to a carrier function, 
for which an analytical Laplace transform exists in order to generate an impedance 

15 spectrum and storing the result of calculation. 

The voltage/current output unit 30 includes first and second filters 310 and 
320 for removing noise by filtering from the voltage and current detected by the 
galvanostat/potentiostat 20; a digital/analog converter 330 for converting the digital 
bias voltage and current output from the control means into an analog bias voltage 

20 and current; first and second subtractors 340 and 350 for respectively subtracting the 
bias voltage and current from the digital/analog converter 330 from the voltage and 
current output from the first and second filters 310 and 320; and first and second 
amplifiers 360 and 370 for amplifying output signals of the first and second 
subtractors 340 and 350 and inputting them to the analog/digital converter 40. 

25 The control means 50 includes CPU 510 for controlling a Laplace transform 

impedance measurement for the electrical or electrochemical device 10; an 
impedance calculation means 520 for measuring and computing the impedance 
spectrum of the electrical/electrochemical device 10 from the output signal of the 
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analog/digital converter 40 under the control of CPU 510; and an input/output device 
530 for inputting the output signal of the analog/digital converter 40 and outputting 
the digital bias voltage and current to the digital/analog converter 330 of the 
voltage/current outputting means 30. 
5 The thus-configured carrier function Laplace transform impedance 

spectrometer is designed to operate the galvanostat/potentiostat 20 which applies an 
excitation signal such as a constant current or constant voltage to the 
electrical/electrochemical device 10 in a current pulse method or voltage pulse 
method, or interrupts the excitation signal applied to the electrical/electrochemical 

10 device 10 in a current interrupt method, or connects a constant load to the 
electrical/electrochemical device in a constant load method and finally detects the 
response signal, namely, the time-variation of the voltage or current. I 
In order to remove the high-frequency component of the voltage and current 3 
response signal detected, the voltage/current output unit 30 performs a low-pass - 

15 filtering and removes bias. That is, the galvanostat/potentiostat 20 performs a low- 
pass filtering the detected voltage and current in the first and second filters 310 and 
320 to thereby remove the high-frequency component. The first and second filters 
310 and 320 outputs signals to the input/output device 530 according to the control of 
CPU 5 1 0 in the control means 50. The signals are converted into the analog signals in 

20 the digital/analog converter 330. The voltage and current bias are individually 
subtracted, and the resulting signals are amplified in the first and second amplifiers 
360 and 370. 

The voltage and current output from the voltage/current output unit 30 are 
converted into digital signals in the analog/digital converter 40 in a desired sample * 
25 extraction speed. The digital signal is input to the input/output device 530 by the 2 
control means 50. The excitation signal and response signals digital-converted in the 
impedance measuring means 520 are fitted to a carrier function which is analytically 
Laplace-transformable by and thus the impedance function of the frequency region is 

7 
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obtained by using the parameter values obtained through the fitting. 

FIG. 2 is a signal flow chart showing a method of measuring the Laplace 
transform impedance. 

Hardware is first initialized and then a measurement mode is selected. 
5 If the selected measurement mode is current pulse method or voltage pulse 

method, the constant current or voltage is applied as the excitation signal to the 
electrical/electrochemical device 10. Variation of the voltage V(t) or current I(t) of 
the electrical/electrochemical device 10 resulting from the application of the 
excitation signal is Laplace-transformed into V(s) and I(s) using carrier functions. 
10 The impedance function V(s) and I(s) is calculated from as Z(s)=V(s)/I(s). The error 
function of the calculated impedance Z(s) is also calculated and stored. 

If the selected measurement mode is a current interrupt method, the input 
signal of the electrical/electrochemical device 10 is interrupted after pre-specified 
time. The variation of the voltage V(t) and current I(t) of the 
15 electrical/electrochemical device 10 is Laplace-transformed to V(s) and I(s) using 
carrier function, and serial capacitance C scr is determined from V(t<0) before 
interruption. With the Laplace transformed V(s) as well as I(s) and C scr? Impedance is 
calculated as follows: Z(s) =V(s)/I(s) + l/sC scr . The error function of the impedance 
Z(s) is also calculated and stored. 
20 If the selected measurement mode is a constant load method, switch is turned 

on to thereby connect the constant resistor R with the electrical/electrochemical 
device 10. The variation of the voltage V(t) and current I(t) of the device 10 is input 
and thus the current I(t) is Laplace-transformed into I(s) using carrier function. From 
the Laplace-transformed I(s) and the constant resistance R, the impedance Z(s) is 
25 calculated as follows: Z(s)=Eo/sI(s)-R, where , Eo = V(t=0). The error function of the 
calculated impedance Z(s) is also calculated and stored. 

Hereinafter, the method of and apparatus for measuring carrier function 
Laplace transform impedance will be described in detail. 
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In applying the excitation signal such as current and voltage or load to the 
spectrometer, if the spectrometer is once excited through other step, we have to wait 
for more than the set time to apply the minimal excitation signal to thereby stabilize 
the system for testing. 

5 The rise time of the signal for controlling the constant voltage and current 

excitation signal should be shorter than 1/2 of the sampling time (ts), and the 
sampling time is related to the maximum frequency f max as the following equation 1 . 
[Equation 1] 

•-TT 

J max 

10 

The magnitude of the applied excitation signal should be determined within 
range that the response of the system for testing remains linear. 

For example, in the electrochemical system like battery, the product of the 
applied current and time should be less than 1/10 (generally, 1/20) of the total 
15 capacity of the battery. 

Generally, if the confidence interval is wide enough in the low-frequency 
band, the size of the excitation signal should be reduced, while the response signal 
- should be enlarged enough to obtain adequate signal to noise ratio. 

The applying time of the excitation signal with respect to the minimum 
20 frequency f min is given as . 

./min 

In case the internal properties like the state of charge is changed by the 
application of the excitation signal as in the secondary battery, and thus it should be 
recovered to the original state after the impedance measurement, an inverse excitation 
signal to the desired measurement can be applied. For example, if the excitation 
25 signal of a constant current is used in the impedance measurement, the current of 
opposite sign is applied for same time and thus the state of the battery can be 
recovered. Here, the measuring of the response signal is not required. 

The measurement method of the constant load impedance can be applied to 
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the energy storing devices such as the charged capacitors, and batteries. 

The resistor used as the constant load should be selected considering the 
maximal variation of the measured voltage. 

In this case, the excitation signal for the recovery of the original state can not 
5 be applied, but as the discharging amount generated during the measurement is small, 
it can be used for the testing of the device. 

Moreover, since the apparatus used for exciting the signal does not apply the 
signal exactly the same as the desired excitation signal, we need to measure the 
applied excitation signal. For instance, the potentiostat, a constant voltage controller, 
10 and the galvanostat, a constant current controller, generally have a limited rise time. 
And if the step type pulse signal is excited by using the apparatii, the actually applied 
signal is a mixture of the pulse and ramp rather than the ideal step pulse. 

Therefore, this kind of problem can be corrected by measuring the actually 
applied excitation, but the oscillation at the moment of signal excitation or a long rise 
15 time may affect on the quality of the measurement result in the high frequency region. 

In order to avoid the above problems, we may use the current interruption 
method. 

This method is to interrupt the current after a predetermined time for applying 
the current and thus to determine the response signal. 
20 In determining the response signal, a low-band pass filter is used in order to 

remove an aliasing effect. 

The sampling time is, as shown in equation 1, selected by the maximum 
frequency f max to be measured, and the number of the measuring data points is 
required as many as f max /f min to thereby memorize the measuring data with respect to 
25 the entire impedance spectrum from the minimum frequency f min to the maximum 
frequency f max . 

For instance, in case of extracting the 16 bit data, the required capacity of the 
memory device is 2 Mb for the measurement of 1 mHz to 1 kHz. If such amount of 
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data information is stored in RAM, the cost for fabricating the measurement 
apparatus is thus increased. 

Carrier function Laplace transform impedance measurement allows the 
extraction of data by logarithmic sampling in time interval without influencing the 
5 quality of the measurement. 

In case of minimizing the number of the measuring data points by using the 
logarithmic sampling time interval, the number of data points may be 20-100 in case 
of the frequency range of 1 mHz to 1 kHz. 

In order to reduce the influence from the noise during measurement, generally 
10 a thousand data points are enough and the capacity of the memory device required for 
the digital/analog converter employing the logarithmic sampling measurement 
method is about 2 Kb. 

According to the method of measuring the Laplace transform impedance, the 
carrier function such as following equation 2 is used for fitting of the response signal 
15 in the time interval with respect to the current/voltage pulse and constant load 
excitation signal. 
[Equation 2] 

where the selection of the carrier function is based on a principle of electrical 
20 engineering that the response of the arbitrary linear circuit with respect to the impulse 
excitation time can be expressed as the equation 2, and 1/ r / corresponds to zeros 
of the circuit. 

The arbitrary form of excitation signal can be expressed as sum of a plurality 
of impulses so that this principle can be generalized to the any excitation signal. 
25 To fit a specific signal, it is useful to fix the above zero values by the value 

distributed by log function from 1/ts to 2 f min . 
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The total number of zeros is determined according to the desired precision, 
and for example, during fitting, we may use a method where the number of the zero is 
automatically increased until the size of the least square of the error may be within 
the error tolerance range. 
5 If such constant zeros are used, the values k t of the parameters in the equation 

2 can be calculated by using a linear regression. This calculation can be performed 
faster than other fitting and is advantageous in that it does not use an initial parameter 
conditions. 

The carrier function without zero corresponds to the ideal pulse and as the 
10 zeros are added, the function deviates from the ideal pulse. 

Laplace transform of the carrier function in equation 2 is given by equation 3 
with respect to the complex parameter s. Here, V(s) is a voltage signal carrier 
function; I(s) is a current signal carrier function. 

15 [Equation 3] 

{a)F(s) = !{/(/)} = Xr—^- 

,=o 1 
s + — 

(b)V(s) = L{v(t)} 
= £{/(/)} 

Impedance function is determined from the following equations 4 to 6 out of 
the correlation of the specific excitation signal and the response signal in the time 
20 domain. 

[Equation 4] 
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V(s) 

Z ^ = 77T (current/voltage pulse method) 



[Equation 5] 



[Equation 6] 



I(s) 



~7< \ V (S) 1 

^ U) = -7— + — - — (current interruption method) 



Z (s) = — ~ R (constant load method) 



In equation 5, C scr is the value calculated from the slope of the response signal 
10 over time just before the applied signal is interrupted. 

In equation 6, E Q is open circuit voltage before the measurement of the tested 
device, and R is the value of the constant resistor connected to the tested device. 

In case of the constant load measurement in equation 6, the initial state 
corresponds to the ideal potential pulse E^s in the time domain, it does not need the 
15 information about the excitation signal. The actual measurement apparatus does not 
require a function of applying the excitation signal other than connecting with the 
resistor, and thus is configured as a single input channel input device for measuring 
the current response signal. 

By substituting 2 K i v (f min < f < f max ) for parameter s with respect to the 
20 result of calculating each of the impedance functions with respect to the constant 
voltage, constant current or constant load by using the equations 3 to 6, we can get 
the impedance value of the tested device for each of frequencies. 

By using the correlation between the standard deviation of the parameters and 
the linear regression in the time domain, we can calculate the deviation for each of 
25 calculation results in the frequency region, and thus calculate the statistic error of the 
frequency domain data using the statistic measurement out of the following equation 
7. 



13 
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[Equation 7] 

AZ(s) = D(s) T C D(s) 

Where D(s) is column vector made by differentiating the impedance values of 
5 equations 3 to 6 from s; and C is a covariance matrix calculated by linearly fitting the 
measured test data. 

EMBODIMENT 1 : 

Measurement times respectively in the frequency scanning method using FRA, 
10 Fourier-transform impedance measurement method and carrier function Laplace- 
transform impedance measurement method are compared in table 1 . Condition of the 
measuring frequency is ranged from 1 mHz to 1 kHz. 

Laplace-transform impedance measurement employs the current pulse method 
and its apparatus is as follows. The current pulse defined at the 
15 galvanostat/potentiostat 20 is applied to the spectrometer, electrical/electrochemical 
device 10. The response signal of the spectrometer, namely, the terminal voltage and 
current are extracted as of digital value from the 2 channel 16 bit analog/digital 
converter 40 through the voltage/current output unit 30. 

The extracted response signal is transmitted to the control means 50 and used 
20 for calculating complex impedance according to the Laplace transform algorithm. 

In the FFT method, the current input signal of superposed sine waves at an 
odd frequency of minimum frequency is applied to battery through the 16 bit 
digital/analog converter and galvanostat. Then, the 2 channel 16 bit analog/digital 
converter measures the voltage and current of the spectrometer, and the complex 
25 impedance is calculated by using digital discrete Fourier transform algorithm. 
[TABLE 1] 



Method of 


Frequency 


Fourier transform 


Laplace transform 
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measurement 


scanning method 


method 


method 


Time required for 
measurement 


1.5 Hour 


33 Min. 


17 Min. 



From the above TABLE 1, we can recognize that the Laplace transform 
method is evidently at least more than two times faster than the FFT method and 5.3 
times faster than frequency scanning method. 

5 

EMBODIMENT 2 : 

Dummy cell consisting of resistors and capacitors is charged by IV and 
connected with the constant resistance 10 Q. Variation of voltage for 1000 seconds is 
measured by time function. 

10 The voltage response signal of the battery is transmitted through the low-band 

pass filter and differential amplifier and then sample-extracted in speed of 1 kHz 
through the 16 bit analog/digital converter. Here, the dummy cell is made by 
connecting RC element where 10 & resistance and 10F capacitor are connected in 
series and other RC element where 1 Q resistance and 1 F capacitor are connected in 

15 parallel. 

The complex impedance spectrum is obtained by transforming the voltage 
data obtained according to the Laplace transform algorithm. 

FIG. 3 shows the complex impedance spectrum obtained by the constant load 
method. The spectrum is accurately identical with the impedance spectrum calculated 
20 from the circuit configuration of the dummy cell. 

EMBODIMENT 3 (COMPARATIVE): 

The impedance spectrum obtained by the carrier function Laplace transform 
method is compared to the impedance spectrum by FFT method. 
25 To the dummy cell configured as the circuit shown in FIG. 4a, current pulse is 



15 
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interrupted after applying current of 0.6mA for 10 minutes. Variations of voltage and 
current are measured for 10 minutes in speed of 10 kHz through the filter and 
amplifier, using the 16 bit analog/digital converter. 

. The thus-obtained current/voltage data is obtained as the linear fitting of the 
5 parameter Kj values of the current/voltage in accordance with the equation 3 and the 
impedance function is calculated according to the equation 4. The measured response 
signal is Laplace-transformed using carrier function, thus obtaining the impedance 
spectrum. 

Frequency of the impedance spectrum ranges from 0.1 Hz to 5 kHz. 
10 The measurement of the impedance spectrum by the FFT method is 

performed by the following steps and apparatus. 

The current input signal of superposed sine waves at an odd frequency of 
minimum frequency 0.1 Hz is applied to battery through the digital/analog converter 
and galvanostat. The current and voltage signals measured by the analog/digital 
15 converter is transmitted to the control means and used for calculating the complex 
impedance, using fast Fourier transform algorithm. 

The impedance obtained by FFT ranges from 0.1 Hz to 18 kHz. Measuring 
time is 20 seconds. 

FIG. 4b compares the impedance measuring data. From the drawing, we can 
20 recognize that the two impedance spectra are identical. 

As shown in the graph, the two impedance spectra are completely identical in 
a high frequency band. But, the impedance spectrum measured by FFT which applies 
a perturbation signal made of non-overlapping frequencies provides less information 
in a low frequency region. On the contrary, the impedance spectra from the carrier 
25 function Laplace transform is advantageous in increasing the measurement resolution 
in both the high and low frequency ranges. 

EMBODIMENT 4 : 
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By making a half cell which consists of an active material LiCo0 2 as an 
operational electrode, the impedance is measured by a current interrupting method. 
The half cell is configured as of the sandwich type three electrode method. For 
reference electrode and counter-electrode is used lithium metal and for electrolyte 
solution is used LiPF 6 in EC/DMC (1:1). 

The impedance is measured in discharging 20 % of the entire capacity in a full 
charged condition. The impedance is discharged with current of 0.203mA for 16 
minutes and 40 seconds. Then, current is cut off for 1 000 seconds to measure the 
response signal, the current interrupting method being repeated. Used is sample 
frequency of 500 Hz for this. The response signal of the voltage/current is Laplace 
transformed, and shown in FIG. 5 as a function of the level of discharge. 

EMBODIMENT 5: 

The impedance spectrum of the primary alkaline battery(manufactured by 
15 Duracell Co. Model name: MN1500) is measured by Laplace transform measurement 
method. Used is the current interruption method where current of 28.5 mA is applied 
for two hours to discharge battery and the applied current interrupted thereafter. 

The response signal is measured for one hour at the sample frequency of 100 
Hz, and 1500 response signals of voltage/current are taken at the same intervals with 
20 respect to the log value of the time. 

The thus-processed response signal is transformed according to the carrier 
function Laplace transform algorithm, thereby obtaining the impedance spectrum. 
The response signal is linearly fitted to the carrier function having an analytical 
Laplace transform to thereby calculate the parameters of the carrier function. 
25 The impedance spectrum is then calculated in frequency domain in accordance with 
parameters for carrier functions. . 

The measured impedance frequency ranges from 1.8 mHz to 3 14.2 Hz. 

The measurement error of the frequency domain impedance function is 
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calculated out of the standard deviation of the parameters of the calculated Laplace 
domain impedance function using the co variance matrix of the time domain fit. 

In FIGS. 6b and 6c is shown a frequency dependent function of the 
measurement error in each section of 95 % confidence interval of the real and 
5 imaginary values of the impedance. 

The error becomes larger in the high frequency band but its size does not 
exceed \Q, thus being considerably accurate result. It is the additional advantage 
provided with the present invention, that the carrier function Laplace transform 
measurement method, determines an error of the resulting measurement. 

10 

EMBODIMENT 6: 

The carrier function Laplace transform impedance measurement is performed, 
using the current interruption method where a current pulse is applied to the lithium 
ion battery with norminal capacity of 1300 mAh (manufactured by Matsushita Co.) 

15 fully charged using constant current at 10 hours at a room temperature. 

The discharge current of 130 mA is applied for 16 minutes and 40 seconds to 
the battery, and then the applied current is interrupted, thereby the response signal 
being measured for 1000 seconds. The response signal is sampled at a frequency of 
500 Hz. The data for current and voltage each having 500,000 data points are 

20 extracted to 2,000 points in the equi-logarithmic spacing in the same time interval. 

The step allows to use the 8 Kb memory by extracting the data at a log time 
interval at the analog/digital converter. The response voltage and current signal are 
converted according to the carrier function Laplace transform algorithm to thereby 
obtain impedance spectrum. 

25 FIG. 7 shows the thus-measured impedance spectrum. The frequency of the 

impedance spectrum ranges from 6.5 mHz to 1.57 kHz. 

Industrial Applicability 
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It will be apparent to those skilled in the art that various modifications and 
variations can be made in the present invention without departing from the spirit of 
scope of the invention. 

Thus, it is intended that the present invention cover the modifications and 
5 variations of this invention. 

For example, for the measurement of the excitation signal and response signal 
for obtaining the impedance spectrum can be used the general analog/digital 
converter instead of the log time interval extracting apparatus. 

Moreover, the low band pass filter is not necessarily used for signal 
10 processing. 

In case of using the high sensitivity wide band analog-digital converter, the 
bias signal compensator and signal amplifier may not be used, which is also the case 
if the bias signal is so small or the size of the response signal is so big within the 
scope maintaining the linearity of the system for testing. 
15 The present invention is described in detail, focusing on the cases when 

applying the constant current, constant voltage in which calculation of analytic 
Laplace transform is relatively simple (current pulse/voltage pulse method), when 
measuring the response signal after interrupting the applied signal (current interrupt 
method) and also when applying the excitation signal of constant resistor (constant 
20 load method), but other excitation signal can be applied. 

Any functions available for analytically transform in the Laplace domain can 
be used as the Laplace transform carrier function. 

It is possible to calculate the parameters of the carrier function by using the 
non-linear fitting with variable zeros, instead of a fixed zero having log time interval. 
25 In a specified case, zeros having other time intervals may be more convenient. 

In addition, if the transfer functions of the test device is of interest rather than 
the individual impedance values at the different frequencies, the impedance function 
obtained from the Laplace domain can be used directly. 
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What is claimed is: 

1. A method of the carrier function Laplace transform impedance 
measurement, comprising the steps of: 

(1) detecting a response signal of an object whose impedance will be 
5 measured; 

(2) obtaining parameters of a carrier function by linearly or non-linearly 
fitting the response signal in step (1) to the carrier function allowing an analytical 
Laplace transform; 

( 3 ) calculating an impedance function in the Laplace domain by using an 
10 analytic relation between the carrier function in time domain and its Laplace transfer 

function using parameters of the carrier function obtained from the fitting results in 
step (2); 

( 4 ) calculating the impedance spectrum in a frequency domain by using the 
parameters of the Laplace domain impedance function calculated in step (3), and 
15 calculating the measurement error of the frequency band impedance function out of 
the standard deviations of the parameters and its correlation; and 

(5) displaying the frequency domain impedance spectrum calculated in step 
(4) and storing the result. 

20 2. The method as defined in claim 1, wherein the above step (1) comprises 

the steps of: 

applying a current/voltage pulse to the object to be measured; and 
detecting the voltage and current based on the voltage pulse and current pulse 
applied as a response signal. 

25 

3. The method as defined in claim 1, wherein the step (1) comprises the 
steps of: 

applying the pulse voltage or current pulse to the measurement object 
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followed by interruption; and 

determining the voltage and current after the interruption as the response 

signal. 

5 4. The method as defined in claim 2 or 3, wherein the measurement 

object is electrical circuit, linear device, non-linear device, capacitor, supercapacitor, 
primary and secondary batteries or fuel cell. 

5. The method as defined in claim 1, wherein in the step (1) a constant 
10 resistor is connected to the measurement object to thereby determine the voltage and 

current as response signal. 

6 . The method as defined in claim 5, wherein the object for measurement is 
capacitor, supercapacitor, primary and secondary batteries or fuel cell. 

15 

7 . A carrier function Laplace transform impedance spectrometer, 
comprising: 

a galvanostat/potentiostat for applying constant voltage, constant current or 
constant load to the electrical/electrochemical device by employing constant voltage 
20 source, constant current source or constant load element such as a resistor and thereby 
detecting voltage and current of the electrical/electrochemical device resulting from 
the constant current or voltage application; 

a voltage/current output unit for removing noise out of voltage and current 
detected in the galvanostat/potentiostat, also removing bias voltage and current and 
25 outputting the resulting signal; 

a two channel analog/digital converter for converting the voltage and current 
output from the voltage/current output unit into digital signal; and 

control means for fitting the digital voltage and current output from the 
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analog/digital converter to carrier function, which is analytically Laplace 
transformable to thereby generate impedance spectrum and storing the result. 

8. The spectrometer as defined in claim 7, wherein the voltage/current 
5 output unit comprises: 

first and second filters for individually filtering the voltage and current 
detected by the galvanostat/potentiostat to thereby remove noise out of them; 

a digital/analog converter for converting the digital bias voltage and current 
output from the control means into analog bias voltage and current; 
10 first and second subtractors for subtracting the bias voltage and current 

output from the digital/analog converter respectively from the voltage and current 
output from the first and second filters; and 

first and second amplifiers for amplifying output signals from the first and 
second subtractors and inputting them to the analog/digital converter. 

15 

9. The spectrometer as defined in claim 7, wherein the control means 
comprises: 

CPU for controlling a carrier function Laplace transform impedance 
measuring operation of the electrical/electrochemical device; 
20 an impedance calculating means for calculating the impedance of the 

electrical/electrochemical device using the output signal of the analog/digital 
converter under the control of CPU; and 

an input/output device for inputting the output signal of the analog/digital 
converter and outputting the digital bias voltage and current to the digital/analog 
25 converter of the voltage/current output unit. 
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FIG. 2 
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FIG. 3 
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FIG. 4a 
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FIG. 5 
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FIG. 6b 
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